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Abstract. In this paper we describe Backlund transformations and hierarchies of exact 
O I solutions for the fourth Painleve equation (PIV) 

^ ■ d^w 1 fdwY 3 , , 2 , 2 , /3 

a^■ di^ = 2^Uj + W + + (1) 

K*" ■ with a, (3 constants. Specifically, a nonlinear superposition principle for PIV, hierarchies 

I of solutions expressible in terms of complementary error or parabolic cylinder functions as 

' well as rational solutions will be derived. Included amongst these hierarchies are solutions 

CN ■ of (1) for which a = ±^?t- and /? = — with n an integer. These particular forms arise 

^ in quantum gravity and also satisfy a discrete analogue of the first Painleve equation. We 

Q.^ ■ also obtain a number of exact solutions of the discrete fourth Painleve equation 

-2znxl + (t/ - 3.5-^ - z^)xl + fi'^ 

[Xn + Zn+ 7)(X„, + Zn - 7) 

O ■ where Zn = n6 and t], 5, /x and 7 are constants, which, in an appropriate limit, reduces to 

PIV (1). A suitable factorisation of (2) facilitates the identification of a number of solutions 
which take the form of ratios of two polynomials in the variable Zn and the limits of these 
^ ' solutions yield rational solutions of (1). 

w. 

1. Introduction 

The six Painleve equations (PI-PVI) were first derived around the turn of the century in an 
investigation by Painleve and his colleagues into which second-order ordinary differential equations 
have the property that the singularities other than poles of any of the solutions are independent of 
the particular solution and so are dependent only upon the equation (cf., [20]); this property is now 
known as the Painleve property. There has been considerable interest in Painleve equations over 
the last few years primarily due to the fact that they arise as reductions of soliton equations solvable 
by inverse scattering as first demonstrated by Ablowitz & Segur [2] . Although first discovered from 
strictly mathematical considerations, the Painleve equations have appeared in various of physical 
applications (cf., [1] and the references therein). The Painleve equations may also be thought of as 
nonlinear analogues of the classical special functions though they are known to be transcendental 
since their solution is not expressible in terms of elementary functions. However rational solutions 
and one-parameter families of solutions of the Painleve equations which can be written in terms 
of special functions are known to exist for particular values of the parameters. For example, there 
exist solutions of PII, PHI and PIV that are expressed in terms of Airy, Bessel and parabolic 
cylinder functions, respectively (cf., [18]). 
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Recently there has been considerable interest in intcgrable mappings and discrete systems, 
including discrete analogues of the Painleve equations. Some of these mapping and discrete 
equations arise in physical applications. For example, a discrete analogue of PI (d-PI) arose in 
the study of the partition function in a two-dimensional model of quantum gravity [9,11,13,19]. 
Subsequently a discrete analogue of PII (d-PII) was derived in [25,27] and later discrete analogues 
of PIII-PV (d-PIII-d-PV) were obtained by Ramani et al. [28] using the method of singularity 
confinement [17]; for further details on the derivation of the discrete Painleve equations see, for 
example, [12,16]. One important result of the investigations is that the form of the discrete Pain- 
leve equations is not unique since there exist several possible discrete analogues of the Painleve 
equations. Kajiwara et al. [21] and Grammaticos et al. [15] have derived exact solutions of d-PII 
and d-PIII in terms of discrete Airy and discrete Bessel functions, respectively, in analogue to the 
aforementioned results for the associated continuous Painleve equations. We further remark that 
Lax pairs and isomonodromic deformation problems are known to exist for d-PI [13], d-PII [25] 
and d-PIII [26]. However, at present, there is no discrete analogue of PVI, nor are there Lax pairs 
for the versions of d-PIV and d-PV derived in [15]. 

In recent work [6,7], we have been concerned with the investigation of Backlund transformations 
and exact solutions for fourth Painleve equation (PIV) 



WW 



" - I {w'f + ^w^ + W + 2{z^ - a)w^ +P, ' = d/dz, (1.1) 



where a and (3 are arbitrary constants, together with an examination of various applications of 
these solutions to several physically motivated nonlinear partial differential equations. In [6] we 
demonstrated how all known exact solutions of (1.1) can be categorised into one of three families; 
in two of these solutions can be determined in terms of the complementary error and parabolic 
cylinder functions whilst the third family consists of solutions which can be expressed as the ratio 
of two polynomials in z. 

In §2 we review the known Backlund transformations for PIV (1.1). In particular it is shown that 
there exist three-term recurrence relations, or nonlinear superposition formulae, for PIV, which can 
be effectively used for the derivation of solution hierarchies of PIV involving algebraic manipulations 
alone. In §3 we briefly discuss hierarchies of exact solutions for PIV including one-parameter families 
of solutions expressible in terms of parabolic cylinder functions, solutions expressible in terms of 
complementary error functions and rational solutions. In §4 we show that one family of solutions for 
PIV, the so-called half-integer hierarchy, generate solutions of the discrete first Painleve equation 
and are relevant to problems in two-dimensional quantum gravity. 

In §5 we discuss some new solutions of the discrete fourth Painleve equation (d-PIV) given by 
the three-point, non-autonomous mapping 

, / , X _ -22:„x3 + (?? - 3(5-^ - zl)xl + fi^ 

with the variable Xn to be found in terms of Zn ^ n5 + (. This mapping is identified as the 
discretised version of PIV. This is observed by taking the limit of (1.2) as (5 — 0, with 7 = 1/S and 
rj and /x finite. This process yields (1.1) with the parameters a and (3 in that equation related to rj 
and /J, according to a = ^rj and /3 = —5^^- Since d-PIV reduces to (1.1) in the appropriate limit 
then it can be expected that exact solutions of (1.2) exist which should tend to known continuous 
solutions in the same limit. Exact solutions of d-PIV additional to those described in §5 have been 
derived very recently and a discussion of these may be found in [4] . 
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2. Backlund transformations for the fourth Painleve equation 

The first Backlund transformation for PIV was derived by Lukashevich [23] who wrote PIV (1.1) 
in terms of a system of Riccati equations. Rather than adopt the method given by Lukashevich, 
we follow Gromak [18] and consider the system 

w' = q + 2ezw + ew'^ + 2£wv, v' = p — 2ezv — ev^ — 2ewv, (2.1a,b) 

where = 1, = —2/3 and p = —1 — ae — If v is eliminated between these equations then 

it is easily shown that w satisfies PIV (1.1). However, if w is eliminated then it transpires that v 
also is a solution of PIV though not with the parameters a and f3 as in (1.1), but rather with ai 
and (3i where ai = \ {2e — 2a + 2>e\/ —2(3) and /3i = —\ (l + ae + \^—2f5) . As equation (2.1a) 
may be rewritten in the form 



w' — 2ezw — £w^ - ^/-2P 



it is apparent that if w{z;a,P) denotes a solution of PIV for parameters a and /? then setting 
£ = ±1 in (2.1-2.2) gives rise to two further solutions of PIV. If the transformations from w{z; a, (3) 
to these new solutions are denoted by and then these Backlund formulae can be written 
as 

~±, , , ^± w\z; a, (3) - w^jz; a, (3) - 2zw{z; a, (3) ^ V-2/j 

W {w{z;a,(3)) : = w{z-a ,(3 ) = — , (2.3a) 

lw\z; a, p) 

5± = i (2-2a±3v/^) , = -i (l + a±iv^^)^ (2.3b,c) 

w±/ / / -± «±x w\z\a,l3) +w'^{z\a,l3) + 2zw{z\a,l3) t V^2^ ,^ . . 

W {w{z;a,P}) ■. = w{z;a ,(3 ) = — , (2.4a) 

2,w\z\ a, j3) 

a± = -i (2 + 2a±3y^) , = -i (l - a ± (2.4b,c) 

Thus given a solution w{z]a,l3) ^ of PIV, the solutions u]{z;a^,(3^) and w{z;a^ , (3^) may 
be obtained using (2.3) and (2.4). Moreover, the transformations (2.3) and (2.4) are effectively 
nonlinear differentiation formulae for solutions of PIV since w{z]di^ , {3^) and w{z]a^ , f3^) are 
expressed in terms of w{z;a^f3) and its derivatives. 

Kitaev [22] derived two further sets of Backhuid transformations for PIV using the associated 
isomonodromy deformation representation, which has a regular singular point at A = and 
an irregular singular point of rank 2 at A = oo (cf., [14]). Kitaev showed that the Backlund 
transformations W and W correspond to the regular point at A = and the transformations 
and W^, defined below, are associated with the irregular point at A = oo. If w{z; a, j3) is a solution 
of PIV, then so also are w\z\a.\f3^) and w'^ {z; , (3^) where 

W^^{w{z-a,P)) ■.=w^^{z-Q^^,l3^^) 

at± = I _ 1^ ^ s^/Z2p, /?t± = - 1 (l - a ± | V^)' , (2.5b,c) 

W^^{w{z;a,P)) ■.=w^^{z;a^^,P^^) 

_ , ..^ 2{\ + a±\V^)w{z-a,l3) 

- w[z, a, p) + ^ _ _ , ^/.oaj 

a*^ = -|-iaT|y^, /3^± = -i(-l-a±i^/^)^ (2.6b,c) 
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These two transformations are valid for all solutions w;(z;q, /3) for which the numerators and 
denominators are non-zero. The transformation W^^ is equivalent to the transformations Ti 
(lower sign) and T2 (upper sign) given by Murata [24]. In [6], is was shown that the Backlund 
transformations W^"^ and W^^ are expressible in terms of W"^ and VF^ 

^t± . . „t± 0t±) = I W+W'^iwiz; a, /?)), if 1 - a - > 0, 

^' ^ ^ \w-W'^{w{z;a,P)), ii 1 - a - ^^/^ < 0, 

^t±(,.^t± 0t±) = lw+W^{w{z;a,P)), if 1 + « + 1 ^=2^ > 0, 
^ ^ \W-W^{w{z;a,l3)), ii 1 + a + ^^/^ < 0. 

Consequently the following Kitaev fractional transformations can be derived 

W^(w{^z-a,^)) 
W^{w{z;a,(3)) 

which are three-term recurrence relations, or nonlinear superposition formulae, for PIV (see [6] for 
further details. These results can be used for the efficient derivation of solution hierarchies of PIV 
by use of algebraic manipulations alone. 

Fokas et al. [14] who, using a Schlesinger transformation formulation related to the isomono- 
dromy method, also deduced a number of Backlund transformations for PIV. However, as for the 
Backlund transformations W^^ and W^^ above, these are expressible in terms of and [6] . 



3. Solution Hierarchies for the fourth Painleve equation 

3.1. Introduction. In this section we shall discuss various one-parameter families of solutions 
and rational solutions for PIV. In order to characterise such families we start by recalling that the 
Riccati system associated with PIV is given by (2.1). Eliminating either v 01 w yields PIV. If u = 
in (2.1) then p = so that q = —2{a£ + 1) and hence 

(3 = -2{l + aef, (3.1) 

which we shall refer to as the one-parameter family condition. One-parameter solutions are found 
by solving the Riccati equation, 

w' - lezw - £w^ - 2(1 + ae) = 0. (3.2) 

Applying the linearizing transformation w = —eu'/uio this form yields the Weber- Hermite equation 

u" - 2ezu' - 2{a + e)u = 0. (3.3) 

Further, if we make the transformation u{z) = r]{^) exp (^e^^), with ^ = V2z, then we obtain the 
parabolic cylinder function equation (cf., [3]) 



(3.4) 
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with u = —a — 1(1 + e). If z/ Z"*", then this equation has the general solution 

where A and B are arbitrary constants and D^{(^) is the parabohc cyhnder function (of., [3]). If 
u G Z"*", then Dn{C) = He„(^)exp where He„(^) is the Hermite polynomial of degree n 

given by 

He„(0 = (-l)"exp {^e) ^ [exp {-'^e)] • (3-5) 

In general, the one-parameter families of solutions for PIV expressible in terms of parabolic 
cylinder functions are given by 

w{z; -(z. + 1), -2^.2) = -V, -2{v + If) = -2z - (3.6a,b) 

w{z; v + l, -2^.2) = ^(^; ^, -2(z. + If) = -2z - (3.6c,d) 

where (f)u{z) satisfies the Webcr-Hcrmite equation (3.3) with e = 1 and a = —v — 1 and tpu{z) 
satisfies (3.3) with e = — 1 and a = + 1 (cf., [10,23]). If = n, with n a positive integer, then 
4>n{z) and tpniz) are polynomials of degree n expressible in terms of the Hermite polynomial IIe„(^) 
(3.5). 

3.2. The Complementary Error Function Hierarchy. If the parameters in (3.2) are a = 1 
and e = —1, then (3.3) may be integrated to yield 

u{z) = B - Aevfciz), (3.7) 

where A and B are arbitrary constants and erfc(z) is the complementary error function given by 

erfc(z) = -= / exp(-t2) dt. (3.8) 

Hence we obtain the exact solution of PIV 

(see also [18]). 

Applying the Backlund transformations W^, W^, W^^, and W^^ to (3.9) yields a hierarchy 
of solutions of PIV of the form 

w{z; mi + 1, -2ml) = R{z, '^{z)), 

where rrii and m2 are integers such that ^(m-i + m2) G Z"*" (i.e., either mi and m2 are either both 
even integers or both odd integers), and R{z,'^{z)) is a rational function of its arguments; see 
Tables 3.2.2 and 3.2.3 in [6] for further details. For example, applying the transformation 1^*+ to 
(3.9) yields 

Wi+ {w{z; 1, 0)) = w{z; -2, -2) = w{z; 1, 0) + ~ ^ ^ , 

W+ {w[z; 1, 0)) 

and so since 

W+ (wiz; 1, 0)) = w{z; 0, -2) = -2z - w{z; 1, 0) = -2z - *(z) 
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then 

w{z;-2,-2)=^{z)- 



2z + 



We remark that the subset of this hierarchy given by w{z; 2n + 1, 0) correspond to the bound 
state solutions discovered in [8] (see also [5] ) , which have the property that they decay exponentially 
as z ^ iboo and are analogues of the bound state solutions for the linear harmonic oscillator (cf., 
[29]). 

3.3. The Half-integer Hierarchy. If the parameters in (3.2) are a = ib^ and /? = — i, then 
under the transformation w{z) = ±9'^{z) / 9±{z) with 0±{z) = r/(^) exp {Tj^'^) and ^ = V2z, it 
follows that ri{^) satisfies (3.4) with u = — |. Hence ry(^) = ^-D_i/2(0 + -S-D_i/2(— 0; "^^^^ A and 
B arbitrary constants, and so we obtain the following solutions of PIV 

w{z-±\-\) = ^ = -z±Q{z), (3.10) 

where 

^/2 \AD-^,2{V^z) - BD-^/2{-V2z)\ 

Qiz) = -z + --^ Ul^ 1 ^/^^ . 3.11 

AD_y2{^z) + BD_^j2{-V2z) 

Applying the Backlund transformations and (2.3, 2.4) to w{z] — i) yields a hierarchy 
of solutions of PIV which have importance in connection with quantum gravity (see §4 below for 
further details). 

3.4. Rational Solution Hierarchies. It is easily verified that w = l/z satisfies PIV (1.1) 
with the parameter choices a. = 2, (3 = —2 and this is a simple example of a rational solution of 
PIV. Two families of rational solutions for PIV take the forms 

w{z-a,l3) = w{z;a,^) = -2z + ?f-^, (3.12a,b) 

where Pfn{z) and Qm{z) denote some polynomials of degree m consisting of either entirely even or 
else entirely odd powers of z. Murata [24] has shown that (1.1) admits unique rational solutions 
of type (3.12a) or (3.12b) if the parameters are of the form 

{a,p) = {±k,-2{l + 2n + k f); k,neZ, n < -1, k > ~2n, (3.13a) 
{a,P) = {k,-2{l + 2n + kf ); k,neZ, n > 0, k > -n, (3.13b) 

respectively. A third family of rational solutions of PIV is characterised by 

w{z;a,P)=-lz+^^^^, (3.14) 



where {a, f3) = (ni,— 1(1 + 3n2)^), with ni and n2 either both even or both odd integers. An 
extensive discussion of the important properties of these rational solution families together with 
tables containing the first few solutions in each of these three hierarchies is given in [6]. 
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4. Exact solutions of the discrete first Painleve equation 

It is well known that the Korteweg-de Vries (KdV) equation, which possesses a similarity 
reduction to PI, has several discrete versions; one of these is the Kac-Moerbeke equation 

Wn,t = -Wn{Wn+l - W^-l), (4.1) 

which reduces to the KdV equation in the appropriate continuous limit. Fokas et al. [13] 
have demonstrated that in the case of (4.1) a similarity solution is characterised by solving it 
simultaneously with the discrete equation 

n = 2tWn + WniWn+l + Wn + Wn-l) (4.2) 

and the continuous limit Wn = — §^(1 ~ 26^77(2;)), n = — (l + ze^), as e ^ maps (4.2) to 
PI (see [19]). Thus equation (4.2) is a version of discrete PI and it has arisen in the theory of 
two-dimensional quantum gravity (cf., [13,25]). We can use (4.1) to substitute for WnWn+i and 
WnWn-i in (4.2) and these two equations can be rewritten as 

n - Wn,t - "itWn - n + Wn,t - "^tWn - wl 

Wn+1 = , Wn-l = '— . (4.3a,b) 

We observe that (4.3) is essentially the same Riccati system as (2.1) since (4.3a) becomes (2.1a) 
when e = 1 and Wn, Wn+i, n and t are replaced by v, w, p and z respectively and (4.3b) becomes 
(2.1b) if Wn, Wn-i, n and t are identified with w, v, —q and z, respectively. If n is replaced by 
n + 1 in (4.3b) and Wn+i eliminated by using (4.3a) it follows that Wn{t) satisfies 

WnWn,tt = \ {wn,tf + + ^tw\ + 2 (t^ + in) - \r? , (4.4) 

which is precisely PIV with a = —\n and /3 = Thus the solution of d-PI (4.4) can be 

expressed in terms of solutions of PIV with parameters a = — |n, (5 = —^n^ and n G Z. In [6] we 
showed that the first few solutions of (4.4), and hence also (4.1) and (4.2), are 



w±i{t) = wit; ±i, -i) = -tT e(t), (4.5a) 

[@{t)±t][e^{t) -t^±iy 



■2' 2^ 

wMt) = w{t; ±1, -2) = ±_yL_^, (4.5b) 
w±3{t) = w{t;±l-^) = - nr^^A/ .0 , .1 . (4.5c) 



with Q(t) as given in (3.11). We remark that Fokas et al. [13] obtained the solution Wi{t) by solving 
an associated Riemann-Hilbert problem for PIV; the derivation given here is much simpler. 
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5. Rational solutions of the discrete fourth Painleve equation 

Tamizhmani et al. [30] noted that if we set rj = 3(5~^ + 7^ — 2(a^ + 6^) and /x = — 6^, then 
the d-PIV equation (1.2) can be factorised as 

, ^ ^_ {xn + a + b){Xn + a-b){xn-a + b){xn -a-h) 
(x„+i + x„)(a;„_i + x„) - (x„ + z„ + 7)(x„ + z„-7) ' ^ 

and the parameters a = jT] and /3 = — ^/i^ arc now given by 

« = is-' + - i(a^ + ?*'), = -i(a^ - &')'• (5.2) 

Simple rational solutions of (5.1) can be found with x„ proportional to Zn', these are 

Xn = -2zn, a=i(5 + 7, 6=|(5-7, (5.3) 
a;n = -|2„, a=i(5 + 7, 6=i(^-7. (5.4) 

In the hmit as (5 ^ 0, with 7 = 1/S, these discrete solutions tend to w{z;0,—2) = —2z and 
w{z;0, — ^) = —^z, respectively, which are the first members of the PIV hierarchies typified by 
(3.12b) and (3.14). If solutions of (5.1) proportional to l/z^ are sought then it is found that 

^^^_S{S±j)_^ a = |<5±7, 6=i<5±7, (5-5) 

and so, in the limit as S 0, with 7 = 1/5, we have ,x„ w{z;±2,—2) = ±.\/z and these 
continuous solutions are members of the family described by (3.12a). 

More complicated rational solutions of d-PIV (1.2) can be deduced by rewriting (5.1) as the 
pair 

_ {xn + a + h){x^ + a-h) _ (a;„ - a + b)(x„ - a - 6) 

[Xn + ± 7) [Xn + T 7) 

which are discrete analogues of the Riccati equation (3.2). Tamizhmani et al. [30] speculated that 
such a formalism ought to lead to discrete rational solutions though they did not present any such 
solutions. It is a routine calculation to verify that equations (5.6a) and (5.6b) are compatible if 
and only if a = |(5 it 7. [We remark at this stage that (5.6) is not the only possibility for the 
splitting of (5.1). Easy generalisations of (5.6) include the multiplication of the right hand sides 
of (5.6a) and (5.6b) by constants C and 1/C respectively or the taking of the factors in different 
pairings. However, it can be shown that in either of these cases we obtain an incompatible couple 
of equations so that the choice of factorisation (5.6) is not as specialised and restrictive as it might 
first appear.] 

If a = ^(5 lb 7, then we can seek solutions of (1.2) by finding solutions of the simpler form (5.6b) 
(or (5.6a) would do equally well). If we write Xn = Pn/Qn then (5.6b) can be recast as 

Kn+l-Pn+l = (^^n ± 7 + S)Pn - fJ'Qn, (5.7a) 
l^n+lQn+1 = -{Zn T l)Qn - Pn, (5.7b) 

where the 'separation' parameter k^+i can depend both on n and Zn (recall that /j, = a' — b'). 
However, if for the moment we take Kn+i = 1) then eliminating P„ between (5.7a) and (5.7b), 
setting 7 = 1/(5 and taking the limit as (5 — shows that Qn — q{z) where q{z) satisfies 

q" = {z'+^- l)q. (5.8) 
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If in this equation we set q{z) = r){^) with ^ = ■\/2z then we obtain the parabohc cyhnder equation 
(3.4) with u = — ^/U. 

If the variable Qn is eliminated from equations (5.7) and the usual limit taken then it follows 
in a manner similar to that outlined above that rational solutions x„ = Pn/Qn of (1.2) exist which 
tend to a function of the form — He^(^)/Hej„(0 whenever /j, = —2m, with m G N. As illustrated 
in §3.1 above, solutions of PIV (1.1) can be expressed in terms of parabolic cylinder functions 
when the parameters a and /? take certain values. Motivated by these comments concerning the 
(continuous) PIV case, we return to equations (5.7) for the discrete situation. If we write 

(P„, Qn) = (A„, S„) X i-Srr (^ ^n-rnSTl ^ ^ 9^ 

where T{z) denotes the usual Gamma function, and choose = 1, then we obtain the pair 

{zn -mST l)An+i + {zn + S± j)An = /x5„, (5.10a) 
{zn -mST l)Bn+i - {z„ T 7)-Bn = An. (5.10b) 

For m G N we can find exact solutions of these equations with An and Bn taking the forms of 
polynomials in consisting of either only even or only odd powers, and of degrees m — 1 and m 
respectively. The first few solutions in this hierarchy are 

m = l, xn = (5.11a) 

"^-A Xn- 2zl-5{25±^y ^^-"^^ 

?>6{25±^) [2zl-5{2,5±^)\ 
m = 3, Xn = 7^ ..g. , o M ' (5.11c) 

Zn [2zi -6(86 ±37)J 

4,5(5,5 ±27) [24- ,5(11,5 ± 37)] 
' 44 - 45(105 ±37)4 + 3^2(3,5 ±7)(45± 7)- ^ ^ 

5S{36 ± 7) - ± 37)4 + 3<5^(3(5 ± j){A6 ± 7)] . . 

' z„ [44 -205(4,5 ±7)4 +(52(256,52 ±12575 +1572)] ' ^ 

In each case the corresponding value of jj, in (5.10a) is = —6m [{m + 1)5 ± 27] and this leads 
to the respective values of parameters a and (3, as defined by (5.2), given by 

a = 1(5-2 - 7^) T (m + 1)75 - ^5^ [2m(m + 1) + 1] , p = -\m^6^ \{m + 1)5 ± 27]^ . (5.12) 

We emphasise at this stage that the discrete solutions (5.11) are exact and are valid for any 6 and 
7. We can recover continuous solutions by letting 7 = 1/5 and taking the limit as 5 ^ which 
yields a = +(m + 1) and /? = —2rr?. Then from (5.11) we obtain 

1 4z 3(2^2 ic 1) 

u;(z;t2,-2) =T-, ^;(z; +3, -8) = T , ^ w;(z; +4, -18) = +- ^ ^ ' 



z' '-v-'^"' ^2z2^l' "-v-'^-' -I ^z(2z2^3)' 

8^(2^2 ^3) , , 5(4^^ + 12^2^3) 

tt;(2:;+5,-32) A.o 9 , o > ^z; +6, -50) = +- ^ ^ ^ 



4^4^1222 + 3' ^ ' > ^2(4z4 + 20z2 + 15)' 

It is then clear that solutions (5.11) can be thought of as discrete analogues of the (continuous) 
solutions of PIV taking forms given by (3.6a) and (3.6d), with v = m. 
If in (5.7) we let 

(P„, Qn) = {An, Bn) X 5" F ( ' rn6 ± j \ ^ 
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then we obtain 

{zn - mS ±-f){Ar,+l - An) = -fiBn, (5.14a) 
{zn -md± 7)5„+i + {zn =F l)Bn = -A^. (5.14b) 

As previously, exact polynomial solutions of these can be derived for integer m except that now 
An and B„ are of degrees m + 1 and m respectively. Then the first few solutions of this type are 

m = 0, Xn = — 2z„, (5.15a) 

2zl - ()(() ^7) . , X 

m = l, Xn = ^ (5.15b) 

2z„[24-J(5,5t37)] 

"^ = '' = 24-^(2^^7) ' ^ 

m-3 X - 44-4^(7Jt37)4 + 3<^^(3^T7)(2c^T7) 

z„ \2zl - m T 37)] ' ^^^^^^ 

2z„ [44 - 20^(3,5 T 7)4 + (146J^ ^ 957(5 + 157^)] 
' 44 - 4,5(10(5 T 37)4 + 3(52(4,5 T7)(35t 7) ' ^ 

The corresponding values of the parameters a and /? are given by 

a = |((5-2 - ± rn7,5 - |,5^ [2m(m + 1) + 1] , /3 = -^(m + 1)^(5^(771(5 T 27)^ (5.16) 

and in the limit (5 — > with 7 = 1/6, solutions (5.15) reduce to the functions 

1 4z 
k;(^;0,-2) = -2z, w{z; ±1,-8) = -2z ^ -, w{z; ±2, -18) = -2z ^ 



/ , „ ^ 3(2^2 ±1) , , . 8z(2z^ ± 3) 

-»;(z; ±3, -32) = -2z T L o . w;(z; ±4, -50) = -2^ ^ 



z(2z2±3)' V ' ' ^ ^ 4^4 ±12^2 + 3- 

These are all members of the so-called ' — 2z' hierarchy of rational solutions for PIV (1.1) and are 
of the form (3.6b) or (3.6c), with v = m. It is straightforward to obtain further exact rational 
solutions of d-PIV by solving the pairs (5.10) or (5.14) for higher values of m. 

We note here that although we have found some rational solutions of d-PIV (1.2) it can be 
anticipated that further such solutions exist which are not derivable using the procedure described 
above. This deduction follows from the observation that rational solutions of PIV (1.1) are possible 
for all parameter values as described by (3.13) and that so far we have discrete solutions which, 
in the appropriate limit, tend to the continuous solutions of type (3.6). Therefore, there should 
be discrete counterparts to the remaining continuous rational solutions. These can be constructed 
by appealing to some Backlund transformations for d-PIV which were given by Tamizhmani et al. 
[30]. These authors presented a sequence of transformations which, given a solution of (5.1) with 
parameters a, b and 7, yields a further solution of the same equation but now with parameters 
a + S, b and 7. In short, suppose that x„ is a solution of d-PIV with parameters a, b and 7. Then 

~ _ XnXn+l + XnjZn + a) + XnjZn - a) + b'^ - 

Xn — ~ 5 {'^■^') 

Xji -\- Xji-\-l 

with = Zn + also satisfies d-PIV but now with parameters a = a + ^S, 5 = 7 and 7 = 6. 
The subtlety with transformation (5.17) is that Xn is defined on a lattice of points which is offset 
by ^5 from the original. In order to obtain a solution valid at points coinciding with the initial 
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lattice it is therefore necessary to reapply (5.17) whicli lias the overall effect of raising the value of 
a by 5 whilst keeping b and 7 invariant. Clearly M applications of this sequence will increment a 
to o + Md and leaves the other two parameters unchanged so that from a starting solution with 
parameters a = b = (N + ± 7, for TV" = 0, 1, 2, . . . (which are the parameter values in 

(1.2) appropriate to the family of solutions whose first few members are listed in (5.11)) we can 
obtain a solution with a = (M + |)(5 it 7 and b = (A/" + |)(5 ± 7, or, in terms of the parameters a 
and (3 as given by (5.2), 



a = - T ^^(^ + N + l)+[^{M-N)iM + N + l)-{M+ ^f] 5\ (5.18a) 



(3 = -\{M - Nf5^ [{M + N + 1)6 ± 2-ff . (5.18b) 



In the usual limit our discrete rational solutions will then tend to continuous ones with associated 

parameters a = ^{M + A'^ + 1), (3 = — 2(M — A^)^. Now M and N can be chosen so as to force 
these parameters to coincide with the form (3.13a) for any k and n in the permissible ranges; in 
this way we have a mechanism for deducing discrete analogues of all those rational solutions of PIV 
(1.1) which take the form Pn-i{z)/Qn{z). 

A similar argument can be applied to the discrete solutions (5.15) and this demonstrates that 
given these results then suitable application of the Backlund transformations contained in [30] 
will generate another set of exact solutions of (1.2). These comprise the discrete counterpart to 
the '— 2z' rational hierarchy of (1.1) which itself is characterised by the parameter values given in 
(3.13b). 

We remark here that of the three hierarchies of rational solutions for PIV (1.1), the simple 
procedures outlined above have yielded discrete analogues of only two of these; no solutions 
corresponding to the —^z family (3.14) have been found. The reason for this is that although 
the full discrete equation (5.1) admits the solution x„ = —^Zn with a = ^5 =b 7, the splitting of 
this equation into the two Riccati-like forms (5.6) gives a pair of equations whose compatibility 
requires that a = ^^±7. Thus, as noted by Tamizhmani et al. [30], this solution is not linearizable 
through the splitting assumption and thus it is unsurprising that the discrete analogue of the — §2 
hierarchy of solutions cannot be generated in this way. However, we shall now show how use of 
the Backlund transformation (5.17) discussed above can be used to increment the values of a or 6 
by ±(5 and thus lead to some more exact solutions in the discrete — |z hierarchy. (It is noted that 
d-PIV (5.1) is invariant by interchange of parameters a and 6 or by a sign change of either of these. 
Therefore, given that two applications of (5.17) increases a by (5 it is easy to see how a suitable 
changes in sign of a and b or the interchanging of these parameters allows transformations to be 
found that increase or decrease the values of a or 6 by integral multiples of S.) A few examples of 
simpler solutions in this hierarchy are: 

Xn = -^Zn, o=i(5±7, 6 = -i(5±7, (5.19a) 
2^^_5(^T37)^ a = ±7-|5, b = T^+lS, (5.19b) 

x.=-'-^M^^^^, a = ±,-16, b = ^, + l6, (5.19c) 



3(24 ± Sj6 - 26^) ' 

a = ±7 + §(5, 6 = ^7 + U, (5.19d) 



2z„(44- 457^(5^ + ^^) 

■3(24 _3^<5_<52)(2^2 +3^^_ J2)' " ^7^6^^' Tl ' 6>^' 



8z^ - 46{6 ± 37)4 - 2(52(38(52 t G9j6 + 27j^)zl + 96^46^ ^ ll-f6^ + 10-f^6 ^ 



3zn [44 - 45(45 ^ 37)4 ± 97^2(5 ^ 7)] 



a = ±7-|(5, b = T7+^S. (5.19e) 
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In the limit as (5 — 0, with jS = 1, the solutions (5.19) reduce to 



w{z;0,-g) = -^z, w{z;±l,-g) = -^zT-, w{z;±2,--^) = -^z ± 
w{z- 0, - f ) = - 1^ ± w{z-±3,-l) = -h± 



24z g 2 _3(4£^_±4£2j-3)^ 

9'" (2z2-3)(2^2+3)' ±3,-9) - -3Z ± ^^^^^ ^ ^2^2 _ . 



which belong to the — |z hierarchy of rational solutions of PIV (1.1); see Table 4.1.3 of [6] for a 
more extensive list of such solutions. The method of generating the — f-^n discrete rational solutions 
given in (5.19) is intrinsically less satisfying than that employed for the evaluation of solutions 
(5.11) and (5.15) in the other two families; this feature is a direct consequence of the fact that the 
— \zn forms do not arise from a suitable factorisation of d-PIV into a pair of simple compatible 
equations akin to (5.6). Instead the present method is based on a direct implementation of the 
Backlund transformation (5.17), which, for the for more complicated solutions in the hierarchy, 
becomes an increasingly laborious task. However, we believe that the Backlund transformation 
(5.17) will inevitably be needed to make a complete evaluation of each of the hierarchies and 
provide a systematic and efficient procedure for doing this. 
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